EXTENSION OF THE 1/-METRIC 



JOSEPH A. BALL AND AMOL J. SASANE 

Abstract. We extend the i^-metric introduced by Vinnicombe in ro- 
bust control theory for rational plants to the case of infinite-dimensional 
systems/classes of nonrational transfer functions. 



1. Introduction 

The general stabilization problem in control theory is as follows. Suppose 
that i? is a commutative integral domain with identity (thought of as the 
class of stable transfer functions) and let ¥(R) denote the field of fractions 
of R. The stabilization problem is: 

Given P € (F(i?))^^"* (an unstable plant transfer function), 
find C € {¥{R))"^^P (a stabilizing controller transfer function), 
such that (the closed loop transfer function) 

P 
I 



HiP,C) :-- 



(I-CP)-^ [ -C I 



belongs to r{p+^)Ap+^) (is stable). 

Recipes for constructing such C is a central theme in control theory; see for 
example the book by Vidyasagar [24j . 

However, in the robust stabilization problem, one goes a step further. One 
knows that the plant is just an approximation of reality, and so one would 
really like the controller C to not only stabilize the nominal plant Pq, but 
also all sufficiently close plants -P to -Pq- The question of what one means by 
"closeness" of plants thus arises naturally. So one needs a function d defined 
on pairs of stabilizable plants such that 

(1) d is a metric on the set of all stabilizable plants, 

(2) d is amenable to computation, and 

(3) d has "good" properties in the robust stabilization problem. 

Such a desirable metric, was introduced by Glenn Vinnicombe in [25] and is 
called the v-metric. In that paper, essentially R was taken to be the rational 
functions without poles in the closed unit disk or, more generally, the disk 
algebra, and the most important results were that the z^-metric is indeed a 
metric on the set of stabilizable plants, and moreover, it has the following 
nice property in the context of the robust stabilization problem: 
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(P): If the z^-metric between two stabilizable plants Pq and P is less than 
the stability margin fip^^c of Pq ^-nd its stabilizing controller C, then C 
also stabilizes P. 

The problem of what happens when R is some other ring of stable trans- 
fer functions of infinite-dimensional (that is, one time axis and infinite- 
dimensional state space) or multidimensional systems (several "time" axes 
of evolution) was left open. This problem of extending the i/-metric from 
the rational case to transfer function classes of infinite-dimensional systems 
was also mentioned in article by Nicholas Young [26]. In this article, we 
address this issue of extending the z^-metric. 

The starting point for our approach is abstract: we suppose that R is any 
commutative integral domain with identity which is a subset of a Banach 
algebra S satisfying certain assumptions, which we label (Al)-(A4). We then 
define an "abstract" z^-metric in this setup, and show that it does define a 
metric on the class of all stabilizable plants. We also show that it has the 
desired property (P) in the context of robust stabilization for an appropriate 
definition of stability margin /ipo,C- 

Next we give several examples of integral domains R arising as natural 
classes of stable transfer functions of infinite-dimensional and multidimen- 
sional systems which satisfy the abstract assumptions (Al) to (A4). In 
particular, we cover the case of full subalgebras of the disk algebra, the 
causal almost periodic function classes, the class of measures on [0, -|-oo) 
without a singular nonatomic part, and the polydisk algebra. 

The paper is organized as follows: 

(1) In Section [21 we give our general setup and assumptions, and define 
the abstract metric dy. 

(2) In Section m we will show that dy is a metric on the set of stabilizable 
plants. 

(3) In SectionlH we introduce a notion of stability margin /ip^c ^-nd prove 
Theorem 14. 6t this implies in particular that if the z^-metric between 
two stabilizable plants Pq and P is less than the stability margin 
f^Po,c of P and its stabilizing controller C, then C also stabilizes P. 

(4) In Section \5\ we specialize R to concrete rings of stable transfer 
functions of various types, and show that our abstract assumptions 
hold in these particular cases. 

(5) The final Section [6] mentions a loose end which is a direction for 
further work. 

2. General setup and assumptions 

Our setup is the following: 

(Al) R is commutative integral domain with identity. 

(A2) S" is a unital commutative complex semisimple Banach algebra with 

an involution •*, such that R C S. We use inv S to denote the 

invertible elements of S. 
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(A3) There exists a map i : inv 5 — )■ G, where (G, +) is an Abehan group 
with identity denoted by o, and l satisfies 

(11) i(a5) = L{a) + L{b) (a, b £ inv S). 

(12) L{a*) = -t(a) (a E inv 5). 

(13) i is locally constant, that is, l continuous when G is equipped 
with the discrete topology. 

(A4) X £ Rn (inv S) is invertible as an element of R iff l{x) = o. 

A consequence of (13) is the following "homotopic invariance of the index" , 
which we will use in the sequel. 

Proposition 2.1. If H : [0, 1] — )■ inv S is a continuous map, then 

c{Hm = LiHil)). 

Proof. The map h, given by t 1— )• L{H{t)) : [0,1] — )• G is continuous. Here 
[0, 1] is equipped with usual topology from M, while G is equipped with the 
discrete topology, given by the metric 

The image of the connected set [0, 1] under the continuous map h is con- 
nected. But the only connected subsets of G are the singleton sets, since G 
is carrying the discrete topology. Hence l{H{0)) = i{H{l). □ 

We recall the following standard definitions from the factorization ap- 
proach to control theory. 

Definition 2.2. 

The notation F(i?): ¥{R) denotes the field of fractions of R. 

The notation F*: If F G i^^^™, then F* G S""^^ is the matrix with the 
entry in the ith row and jth column given by F*j^, for all 1 < i < p, and all 
1 < j < m. 

Right coprime/normalized coprime factorization: Given a matrix 
P G (F(i?))^^™, a factorization P = ND~^, where N,D are matrices with 
entries from R, is called a right coprime factorization of P if there exist 
matrices X, Y with entries from R such that XN + YD = If moreover 
there holds that N*N + D*D = Im, then the right coprime factorization is 
referred to as a normalized right coprime factorization of P. 

Left coprime/normalized coprime factorization: Similarly, a factor- 
ization P = D~^N, where N, D are matrices with entries from R, is called a 
left coprime factorization of P if there exist matrices X, Y with entries from 
R such that NX + DY = Ip. If moreover there holds that NN* + DD* = Ip, 
then the left coprime factorization is referred to as a normalized left coprime 
factorization of P. We note that the existence of both a left and right nor- 
malized factorization P = ND~^ = D~^N for P leads immediately to a 
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normalized double coprime factorization of P, i.e., one has the identity 



'N* D*' 




'N 


-D* 




I 0" 


-b N 




D 


N* 




/ 



(2.1) 



'N 


-D* 




"AT* 


D*' 




I 


0" 


P 


N* _ 




-b 


N _ 







/ 



Since we are deahng with finite matrices over a commutative ring, (|2.1 
imphes also the identity 



(2.2) 



The notation G, G, K, K: Given P G {¥{R))p^'^ with normahzed right and 
left factorizations P = ND^^ and P = D~^N, respectively, we introduce 
the following matrices with entries from R: 

N 



G 



D 



and G = [ -D N 



In this notation the fact that the left and right coprime factorizations of P 
are normalized translates to 



G*G = I, GG* = 
and the identity (|2.2p assumes the form 

GG* + G*G = I. 



I 



(2.3) 



(2.4) 



Similarly, given G G (F(i2))'"^P with normalized right and left factorizations 
G = NcDf^^ and G = b^^^Nc, respectively, we introduce the following 
matrices with entries from R: 



K 



Dc 
Nc 



and K = [ -Nc Dc ] ■ 



The notation S{R,p,m): We denote by '§>{R,p,m) the set of all elements 
P G (¥{R))P^"^ that posses a normalized right coprime factorization and a 
normalized left coprime factorization. 

Remark 2.3. Given P £ (F(i?))P^™ and G £ (F(i?))'"^P, define the dosed 
loop transfer function 

P 



H{P,G) :-- 



I 



(/ - CP)-i [-G 7 ] G {¥{R)) 



(p+m) X (p+m) 



§{R,p,m) C \ P£ {¥{R))P 



(2.5) 



It can be shown (see for example [231 Chapter 8]) that if P G S{R,p,m), 
then P is a stabilizable plant, that is, 

3G G (F(P))™^P such that 
H{P,G) G i?{p+™)x{p+™) 

It was shown by A. Quadrat [18^ Theorem 6.3] that if the Banach alge- 
bra P is a projective- free ring, then every stabilizable plant admits a right 
coprime factorization and a left coprime factorization, that is, the reverse 
containment D and hence equality holds in (|2.5p . 
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We will need a couple of straightforward results on coprime factorizations, 
which we have listed below. The first lemma says that coprime factorizations 
are unique up to invertibles. 

Lemma 2.4. Let P G (F(i?))P^'". 

(1) If P has right coprime factorizations P = NiD^^ = iV2D^^, then 
there exist V,A£ i2™x"^ such that VA = AV = Im, Ni = N2V and 
Di = D2V. 

(2) // P has left coprime factorizations P = D^^Ni = D2^N2, then 
there exist V,A £ M'^'p such that VA = AV = Ip, Ni = VN2 and 
Di = VD2. 

In the case of normalized coprime factorizations, the invertibles can be 
chosen to be unitary. 

Lemma 2.5. Let P G (F(i?))P^™. 

(1) IfP has normalized right coprime factorizations P = NiDY^= N2D2] 
then there exists a U £ Jimxm^ which is invertible as an element of 
R">^x"\ and such that U*U = UU* = I^,Ni = N2U and Di = D2Uj 

(2) If P has normalized left coprime factorizations P = D^^Ni = D2^N2, 
then there exists a U £ RP^p which is invertible as an element of 
RP^'P, and such that U*U = UU* = Ip, Ni = UN2 and Di = UD2. 

Lemma 2.6. Suppose that F G i?""^™-, detF G inv S and i(detF) = o. 
Then F is invertible as an element of R^^"^. 

Proof. Since det F G inv S and i(det F) = o, it follows from (A4) that det F 
is invertible as an element of R. The result then follows from Cramer's 
rule. □ 

We now define the metric d^ on §{R,p, m). But first we specify the norm 
we use for matrices with entries from S. 

Definition 2.7 (|| • ||oo)- Let SOT denote the maximal ideal space of the 
Banach algebra 5. For a matrix M G SP^"^, we set 

||M||oo = max |M(99)|. (2.6) 

Here M denotes the entry-wise Gelfand transform of M, and | • | denotes 
the induced operator norm from C" to C^. For the sake of concreteness, we 
fix the standard Euclidean norms on the vector spaces C™ to C^. 

The maximum in (12. 6p exists since 9Jt is a compact space when it is 
equipped with Gelfand topology, that is, the weak-* topology induced from 
C{S] C). Since we have assumed S to be semisimple, the Gelfand transform 

^ : S ^ S (C C{Tl,C)) 

is an isomorphism, li M £ S^^^ = S , then we note that there are two norms 
available for M: the one as we have defined above, namely ||M||oo, and the 
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norm || • || of M as an element of the Banach algebra S. But throughout this 
article, we will use the norm given by (12. 6p . 



Definition 2.8 (Abstract z^-metric d^). For Pi,P2 G S{R,p,m), with the 
normalized left/right coprime factorizations 

Pi = NiD^^ = Dy'Ni, 

P2 = N2D^^ = D2^N2, 

we define 

dJPi P2) ■= I 11^2^1 lloo if det(G^G2) G inv S and ;,(det(GtG2)) = o, 
' \ 1 otherwise. 

(2.7) 

Normalized coprime factorizations are not unique for a given plant in 
E{R,p,m). But we have the following: 



Lemma 2.9. d,^ given by <\2.7\} is well-defined. 

Proof. This follows from Lemma 12.51 □ 
Lemma 2.10. d^, given by (j2.7p is bounded above by 1. 

Proof. We have IIG2G1II00 <J|G2||oo||G'i||oo- As GJGi = and G2G^ = Ip, 
we see that ||Gi||oo = 1 = ||G2||oo- D 

In Section [31 we will also prove the following. 
Theorem 2.11. d^, given by (|2.7p is a metric on S(-R,p, m). 

We recall the definition of singular values of a square matrix, and a few 
properties which will be needed in the sequel. 

Definition 2.12. If M G C''^'^, then the set of eigenvalues of MM* and 
M*M are equal and the eigenvalues are real. The square roots of these 
eigenvalues are called the singular values of M, and the largest of these is 
denoted by a(M), while the smallest of these is denoted by a{M). 

Proposition 2.13. The following hold for P,Q £ C'^^*^. 

(51) \\P\\=a{P). 

(52) IfP is invertible, then a{P) > 0, and \\P~'^\\ = (o:(P))"^ 

(53) \aiP + Q) - aiP)\ <^(Q). 

(54) a{PQ)<a{P)-a{Q). 

(55) a{PQ)>a{P)-a{Q). 

(56) a{PQ) = a{{P*P)'2Q) = a{P{QQ*)'2). 

Proof. (SI), (S2) follow from the spectral theorem. (S3), (S4), (S5) are given 
in [H Proposition 9.6.8, Corollary 9.6.6]. (S6) can be verified directly using 
the definition of a. □ 

Lemma 2.14. Suppose that A,B e Cp^"* and that A* A + B*B = L Then 
{a{A)r + ia{B)r = 1. 
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Proof. This follows from the spectral theorem. Indeed, A* A = I — B*B, 
and so for all x E C" with unit norm, we have 

{A*Ax,x) = {x,x) - {B*Bx,x) = 1 - {B*Bx,x). 

Thus (-a(A)f = maxMMx, x) = 1 - min {B*Bx, x) = I - {a{B)f. □ 

lkll=i Il2^ll=i 

In particular, we have the following consequence as an application of this 
lemma. (In this article, we often suppress the argument of the Gelfand 
transforms of matrices with S-entries.) 

Lemma 2.15. If Pi, P2 e S{R,p,m), then (ct(G^Gi)) V (^(GaGi))^ = 1 

pointwise on Tl. 

Proof Observing from ([231) and that Gld = I and G2G^+G^G^ = /, 
we obtain ^ ^ 

]^ Vjr2'jr2»-Jri -\- \or-^*jr2VJr2'-J'l — ^ 

pointwise on DJt. An application of Lemma 12.141 now yields the result. □ 

3. di, IS A METRIC 

In this section, we will prove Theorem 12. 11 [ 
Proof (of Theorem 12. lip . 

3.1. Positivity. If Pi,P2 G S{R,p,m), then clearly d^{Pi,P2) > 0. Also, if 
d„{Pi,P2) = 0, then IIG2G1II00 = 0, and so G2G1 = 0. But 

G2Gi = D2{P2 - Pi)Di. 

Thus Pi= P2. Finally, for P G S{R,p,m), it is clear that d^{P,P) = 0. 

3.2. Symmetry. Let Pi,P2 G S(i?,p,m). Since GIG2 = {G^Gi)*, it fol- 
lows that det(G^G2) is invertible as an element of S iff det(G2Gi) is in- 
vertible as an element of S. Using (12), we see that t(det(G^G2)) = o iff 
t(det(G^Gi)) = o. Hence d^{Pi,P2) = d^{P2,Pi). 

3.3. The triangle inequality. Suppose that Pi,P2,Po G S{R,p,m). We 
want to show that dy{Pi, P2) < dy^Pi^Po) + dy{PQ, P2). Since di, is bounded 
above by 1, this inequality is trivially satisfied if either cij/(Pi,Po) = 1 
or d^{Po,P2) = 1. So in the rest of this subsection, we will assume that 
di,{Pi,Po) < 1 and du{Po,P2) < 1- This means that 

(1) det(GtGo) is invertible in S and i{det{GlGo)) = o. 

(2) det{G*oG2) is invertible in S and i(det(GSG2)) = o. 
We will consider separately the following two possible cases: 

1° det{GlG2) G inv 5 and i(det(GtG2)) = o- Then d^{Pi,P2) = ||G2Gi||oo. 
2° -^[det{GlG2) G inv S and i(det(G];G2)) = o]. Then d^{Pi,P2) = 1. 
First, using the fact (|2.4p that GqGq + GqGq = I, we obtain that 

1<J2 — Lri'j'0<-f0^2 + (-riLro(-rOLr2- W'-'-J 
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1° Suppose that det(G^G2) G inv S and i(det(GJG2)) = o. In this case, 
du{Pi,P2) = ||G2Gi||oo- Using (S3) from Proposition 12.131 with 

Q ■= GiGqGoG2. 

and ([33]), we have o:(G^ GoG^Ga) -^(G^Ga) < ^(G^G^GoGa) pointwise 
on dJl. Furthermore, using (S4) and (S5), and rearranging, we obtain 

a(G^G2) > a(G*iGo) • ^(G^Ga) - a(GoGi) • ^(GoGa) (3.2) 

pointwise on 971. Since d,^{Pi,Po) and diy{Po,P2) are both in [0,1], we can 
find a,f3 (which are maps from 9Jt to [0, ^]) such that 

sin a = a(GoGi), 
sin/3 = a(GoG2), 
pointwise on dJt. Then using Lemma 12.151 it follows from (|3.2p that 

^(GiGa) > (cos a) • (cos (3) — (sin a) • (sin/3) = cos(a + /3) (3-3) 

pointwise on Tl. Similarly, define 7 : 9Jt — )■ [0, |] by ^(GiGa) = sin7, 
then ^(G^Ga) = cos 7 pointwise on 9Jl. The inequality ()3.3p now says that 
cos 7 > cos(q + /3) pointwise on Tl. Hence 

sin 7 < sin(Q + /3) = (sin a) • (cos/3) + (sin/3) • (cos a) 
< (sin a) ■ 1 + (sin /3) • 1, 

that is, ^(GiGa) < ct(GoGi) + ^(GoGa) < d^(Pi,Po) + d^iPo,P2) point- 
wise on 9JT. Consequently, du{Pi,P2) = IIG1G2II00 < d,y{Pi,Po) + d,y{Po,P2)- 

2° ^[det(GIG2) G inv S and i{det{GlG2)) = o]. In this case d^{Pi,P2) = 1. 
Let 

A := GlGoG*oG2, and 

r) /~Hf ri* /~i ri 

D .— Lr^LrQLroL,2- 

Using the fact that GJGq and GqG2 are invertible in S"^^"^, it follows also 
that A is invertible in 5"*x™. 

Suppose that ||^~^i?||oo < 1- Then it follows from (j3.ip that 

GIG2 =A + B = A{I + A^^B) 

and so GIG2 is also invertible in 5""^"^. Consider the map H : [0, 1] — t- inv S, 
given by Hit) = det(A(/ + tA-^B)), t G [0, 1]. By Proposition [231 

o = o + o = i{GlGo) + L{G*oG2) = L{detA) = L{H{0)) 
= i{H{l)) = i{det{GlG2)). 

But then we have that det(G*G2) G inv S and L{det{GlG2)) = o, which is 
a contradiction. 
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So our assumption that ||yl~^-B||oo < 1 cannot be true. From the com- 
pactness of 9Jt and the definition of the norm on C"*^"*, it follows that there 
is a G 931 such that a{{A~^'B){(p)) > 1. But then we have that 

1 < a((A-iB)M) < a((AM)-i) . a{B{^)), 

and so 

a{A{^))<a{B{^)). (3.4) 

Thus 

(1 - (a((GoGi)M))2) . (1 - (a((GoG2)M))2) 
= {a{iG{Go){ip)f ■ {a{iG*oG2){ip)? by LemmaMB 

< a{{GlGoG*oG2){ip)f by (S5) in Proposition EH 

< a((GjGSGoG2M)ny i3D 

< (a((GoGi)((/.)))2 . {a{{GoG2){^))f by (S4) in Proposition m 
With 

X := a((GoGi)((^)), and 

y := a((GoG2)M), 
the above says that (1 — x^) • (1 — y^) < x^y^, and so 1 < + y^. Thus 
1 < (a((GoGi)(v9)))2 + (a((GoG2)M))' < K(Po, A))' + K(Po,^'2))'. 
Consequently, 

{d,iPo,Pl) + duiPo,P2)f > {du{Po,Pl)f + idu{Po,P2)f > 1 =idu{Pl,P2)f. 

Taking square roots, we obtain the desired conclusion. 

This completes the proof of the triangle inequality, and also the proof of 
Theorem Elll □ 

4. Robust stability theorem 
In this section we prove Theorem 14.61 

Definition 4.1. Given P G (F(i2))P^'^ and C £ (F(i?))™^P, we define the 
stability margin of the pair {P, C) by 

_ r \\H{P, C7)||-i if P is stabilized by C, 
I^P,c "j^ Q otherwise. 

The number fip^c can be interpreted as a measure of the performance of 
the closed loop system comprising P and C: larger values of /ip^c correspond 
to better performance, with fip^c > if C stabilizes P. 

Proposition 4.2. If P is stabilized by C, then fipc = inf aCK(ip)G(ip)). 
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Proof. We now write P = NM ^ = M for a normalized left /right 
coprime factorization of P and C = NcM^^ = M^^N^ for a normalized 
left /right coprime factorization of C and we set 



G 



N 
M 



G 



-M N 



K 



Mr 



K 



Then we have 
H{P,C) 



P 
I ^ 

NM-^ 
I 



N 
M 



{i-cpy^ [ -c I ] 



Also, G*G 



= G{KGy^K 

[ N* M* 



N 
M 



Nr Mr 



N*N + M*M = I. Similarly, 



KK* 



-N, M, 



-N* 

m: 



NcN* + McM* = I. 



Using (S6) of Proposition 12.131 we obtain for each 99 G 9K that 

a(GM(K((^)GM)-^K(v.)) = a((K(^)GM)-i: 
Thus 



1 



1 



a(K((^)G((^)) 
1 



supa(GM(KMG(v9))-iK(v9)) = sup — 
A^p,c <pem v'got a;(K((^)G((/?)) 

and so ^p,c = inf a;(K((^)G((/?)). 



□ 



Remark 4.3. It is useful to note that 

^P,C < 1 (4.1) 

for any P and C as above. One way to see this is to note that H{P, C) 
is idempotent H{P,C) ■ H{P,C) = H{P,C); this forces \\H{P,C)\\oc > 1- 
Another way to see (j4.ip is to use the formula for f-tp^c in Proposition 14.2 
as follows. Since G*G = I and KK* = i, it follows that a(G) = 1 and 
^{K) = 1. Then it follows from various of the properties singular values 
listed in Proposition 12. 131 that 

a(K((^)G(<^)) < a(K((^)G((^)) < a(K(<^)) • a(G(^)) = 1. 

Proposition 4.4. The following are equivalent: 

(1) C stabilizes P. 

(2) det(K(99)G(v3)) / for all ip £ m and L{det{KG)) = o. 
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Proof. Suppose that C stabilizes P. Then from the calculation done above 
in the proof of Proposition 14.21 we have 

H{P,C) = G{KG)-^K. (4.2) 

But we know that G is left invertible and K is right invertible as matrices 
with entries from R. So from the above, we see that KG £ j^mrxm jg 
vertible as an element of ij*"^™. In particular det{KG) is invertible as an 
element of R and so det(K{ip)G{(p)) / for aU ip £ Tl. Also from (A4), it 
follows that L{det{KG)) = o. 

Nowjuppose that det{K{ip)G{ip)) / for all £ Tl and i{det{KG)) = o. 
Then KG G i?ninv S. From (A4), we obtain that det(i^G) is invertible as 
an element of R, and so we see from (j4.2p that H{P, C) has entries from R. 
So P is stabilized by C. □ 

Proposition 4.5. ^p^c = fJ'C,P- 

Proof. It is not hard to see that G stabilizes P iff P stabilizes C. We have 

KG*GK* + KGG*K* = /, 

GK*KG* + GKK*G* = / 

pointwise on 9Jt. So it follows from Lemma 12.141 that 

(a(KG))2 = 1 - (a{kG*)f = {a{GK)f. (4.3) 

This completes the proof. □ 

Theorem 4.6. IfPo,Pi G S{P,p,m) and C G S{R,m,p), then 

sin~^ HP^^C > sin"^ hp^^c - sin"^(d^(Po, -Pi))- 

Proof. If d,y{Po,Pi) > /^Po,C) then sm~^{d^{Po, Pi)) > sm~ ^ fip^^^c and so 
sin~-^ /iPg^C" — sin~^((ii^(Poi Pi)) < 0. The claimed inequality in the statement 
of the theorem now follows trivially since iJip^^c ^ 0. 

We therefore assume in the rest of the proof that du{Po,Pi) < fJ'Po,c- As 
noted in Remarli [43l UPg^c < 1: hence we must have d^{PQ,Pi) < 1. Also 
f^Po,c = implies that d^{PQ, Pi) < 0, a contradiction to the fact that d^, is 
a metric. Hence fJ-p^^c > 0, that is, C stabilizes Pq. Now 

d,(Po,Pi) = supa((GoGi)(^)) < inf a{{KGo){p)) = m,c, 

and so pointwise on Tl, there holds that a(GoGi) < ct(KGo). But for 
numbers a,b £ (0, 1), 

a < 6 iff K < 



1 - a2 1 - 62 ' 
and so we have 

(a(GoGi))2 ^ (£(KGo)) 



2 



1 - (a(GoGi))2 1 - (a(KGo))2 
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Using Lemma 12.151 and (14.31) , we obtain '^^^"^"'"^ < . Thus 

a(KG^GoGi) <£(KGoGSGi). (4.4) 

But 

KGi = KGoG*qGi + KGIGqGi. (4.5) 

Let yl := KGoGlGi, and S := KGlG^Gi. Using the fact that KGq and 
GgGi are invertible in it follows also that A is invertible in S^^"^. 

Also, from (fO|) . it follows that ||^~^S||oo < 1- Then it follows from (|i3]) 
that KGi = A + B = A{I + A-^B) and so KGi is also invertible in 5"^^™. 
Consider the map H : [0, 1] inv 5, defined by H{t) = det{A{I + tA-^B)), 
t £ [0, 1]. By Proposition EH it follows that H{0) = H{1), that is, 

L{detiKGi)) = L{detiKGoG*oGi)) = i(det(KGo))+t(det(GSGi)) = o+o = o. 

But det(KGi) G R. By (A4) it fohows that det(KGi) is invertible as an 
element of R. Consequently C stabilizes Pi and 

f,P,,c= inf a((KGi)((^)). 

From (j4.5p . we have 

£(kGi) = a(KGoGSGi+KGSGoGi) 

> a(kGoGSGi)-a(KGSGoGi) 

> £(kGo)a(GSGi) - a{KG*o)a{GoGi) 
= sin(sin-^a(kGo) -sin~^a(GoGi)). 

Since sin"-*^ : [—1, 1] — >■ [— §, is an increasing function, it now follows that 

sin~^o:(kGi) > sin~^a(kGo) - sin"^ ct(GoGi). 

Consequently, sin"-*^ /^Pi,c ^ sin~^ fJ'Po,c — sm~^{di^{Po, Pi)). □ 

Corollary 4.7. If Pq,P G §{R,p,m), then 

f^p,c > f^Po,c - du{Po,P). 

Proof. For x,y, z G [0, 1], if sin""*^ x < sin~^ y + sin~^ z. By taking the cosine 
of both sides and using that the cos is a decreasing function on [0, we 

then get — x"^ > y^l — y^Vl — z"^ — yz, which in turn implies that 



{VT^ + yz)^>{l-y^){l-z'). 

Hence < + + 2yzVl — x"^ < + + 2yz -1 = {y + z)"^, which gives 
finally that x < y + z. The claimed inequality now follows immediately from 
the inequality in Theorem 14.61 upon setting x = hpq^c, U = du{Po,P) and 
z = lJ-p,c- □ 
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The above result says that if the controller C performs sufficiently well 
with the nominal plant Pq) and the distance (ijy(Po) P) between the plant P 
and Pq is sufficiently small, then C is guaranteed to achieve a certain level 
of performance with the plant P. So if P and Pq represent alternate models 
of the system (one which is "true" and one which is our nominal model) and 
if d^{Po,P) is small, then the difference between P and Pq can be ignored 
for the purposes of designing a stabilizing controller. 

Another way of stating the result in Theorem 14.61 is that if C stabilizes 
Pq with a stability margin fip^c > and P is another plant which is close 
to Pq in the sense that du{P, Pq) < m, then C is also guaranteed to stabilize 
P. Furthermore, if C satisfies the stronger condition /ip^c > M > m for a 
number M, then C is also guaranteed to stabilize P with a stability margin 
fip^c which satisfies /ip,c ^ sin~^ M — sin^^ m. 

5. Applications 

Now we specialize R to several classes of stable transfer functions and obtain 
various extensions of the i/-metric. Some of the verifications of the properties 
(Al)-(A4) are similar to the section on examples from [22|. 

5.1. The disk algebra. Let 

B := {z G C : |z| < 1}, B := {z e C : \z\ < 1}, T := {z e C : \z\ = I}. 

The disk algebra A(E)) is the set of all functions / : D — t- C such that / is 
holomorphic in ID and continuous on D. Let C(T) denote the set of complex- 
valued continuous functions on the unit circle T. For each / S inv C(T), we 
can define the winding number w(/) G Z of / as follows: 

w(/) = ^(e(2^)-e(o)), 

where : [0, 2-k] — ?• M is a continuous function such that 
/(e**) = |/(e^*)|e^®W, t G [OM- 

The existence of such a Q can be proved; see \23\ Lemma 4.6]. Also, it can 
be checked that w is well-defined and integer-valued. Geometrically, w(/) is 
the number of times the curve t 1— )• /(e**) : [0, 27r] — >• C winds around the 
origin in a counterclockwise direction. 
Recall the definition of a full subring. 

Definition 5.1. Let Ri,R2 be commutative unital rings, and let Ri be a 

subring of R2. Then i?i is said to be a full subring of R2 if for every x £ Ri 
such that X is invertible in R2, it holds that x is invertible in Ri. 



14 



JOSEPH A. BALL AND AMOL J. SASANE 



Lemma 5.2. Let 

R = a unital full subring of ^(D), 
S := C(T), 
G := Z, 
i := w. 

Then (A1)-(A4) are satisfied. 

Proof. (Al) is clear. The involution •* in (A2) is defined by f*{z) = f{z) {z G 
T) for / G C(T). (A3)(I1) and (A3)(I2) are evident from the definition of w. 
Also, the map w : inv C(T) — t- Z is locally constant (that is, it is continuous 
when Z is equipped with the discrete topology and C(T) is equipped with 
the usual sup-norm); see [231 Lemma 4.6. (ii)]. So (A3) (13) holds as well. 
Finally, we will show below that (A4) holds. 

Suppose that f £ RO (inv C(T)) is invertible as an element of R. Then 
obviously / is also invertible as an element of A(D). Hence it has no zeros or 
poles in B. For r G (0, 1), define fr G A{B) by fr{z) = f{rz) {z G B). Then 
fr also has no zeros or poles in B, and has a holomorphic extension across 
T. From the Argument Principle (applied to /r), it follows that w(/r) = 0. 
But Wfr - /lloo — )• as r 1. Hence w(/) = lim w(/r) = lim = 0. 

r—i-l r— >-l 

Suppose, conversely, that f £ RCi (inv C(T)) is such that w(/) = 0. For 
all r G (0, 1) sufficiently close to 1, we have that fr G inv C(T). Also, by the 
local constancy of w, for r sufficiently close to 1, w(/r) = w(/) = 0. By the 
Argument principle, it then follows that has no zeros in B. Equivalently, 
/ has no zeros in rB. But letting r 1, we see that / has no zeros in B. 
Moreover, / has no zeros on T either, since / G inv C(T). Thus / has no 
zeros in B. Consequently, we conclude that / is invertible as an element of 
A(B). (Indeed, / is invertible as an element of C(B), and it is also then 
clear that this inverse is holomorphic in B.) Finally, since ii is a full subring 
of ^(B), we can conclude that / is invertible also as an element of R. □ 

Besides ^(B) itself, some other examples of such R are: 

(1) RH°°(D), the set of all rational functions without poles in B. 

(2) The Wiener algebra P^+(B) of all functions / G A{0) that have an 
absolutely convergent Taylor series about the origin: 

OO OO 

J2\fn\< +00, where f{z) = ^ {z G B). 

n=0 n=0 

(3) 9-"iJ°°(B), the set of / : B ^ C such that /, /(i), . . . , /(") 
belong to H°°{0). Here ff°°(B) denotes the Hardy algebra of all 
bounded and holomorphic functions on B. 

In the definition of the i/-metric given in Definition 12.81 corresponding to 
Lemma 15.21 the || • ||oo now means the following: if F G (C(T))^^"*, then 

||F||oo = max|F(z)|. 
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This follows from ()2.6p . since the maximal ideal space 9Jt of S = C(T) 
can be identified with the unit circle as a topological space; see [20], Exam- 
ple 11.13.(a)]. 

Remark 5.3. RH°°(I])) is a projective free ring since it is a Bezout do- 
main. Also ^(D), M^^(O), or d~'^H°°(D) are projective free rings, since 
their maximal ideal space is D, which is contractible; see |2]. Thus if R is 
one of RH°°{B), A{B), W^+(B) or 9-"iJ°°(D), then the set S{R,p,m) of 
plants possessing a left and a right coprime factorization coincides with the 
class of plants that are stabilizable by |18l Theorem 6.3]. 

5.2. Almost periodic functions. The algebra AP of complex valued (uni- 
formly) almost periodic functions is the smallest closed subalgebra of L°° (R) 
that contains all the functions e\ := e*^^. Here the parameter A belongs to 
M. For any / G AP, its Bohr-Fourier series is defined by the formal sum 

J]/Ae*^^ yGM, (5.1) 

A 

where 

fx := lim / e-'^yf{y)dy, A G R, 

and the sum in (15. ip is taken over the set <t(/) := {A G R | /a 7^ 0}, 
called the Bohr- Fourier spectrum of /. The Bohr-Fourier spectrum of every 
/ G AP is at most a countable set. 

The almost periodic Wiener algebra APW is defined as the set of all AP 
such that the Bohr-Fourier series (|5.ip of / converges absolutely. The almost 
periodic Wiener algebra is a Banach algebra with pointwise operations and 
the norm ||/|| := |/a|. Set 



AP+ = {f GAP\aif) C [0,00)} 

APW+ = {f e APW \a{f )C [0,00)}. 

Then AP~^ (respectively APW~^) is a Banach subalgebra of AP (respectively 
APW). For each / G inv AP, we can define the average winding number 
w{f) G R of / as follows: 

w{f)=nm -l('arg(/(r))-arg(/(-r)) 

T— i>oo 21 \ 

See [m Theorem 1, p. 167]. 
Lemma 5.4. Let 



R 
S 
G 

w. 

Then (A1)-(A4) are satisfied. 



a unital full subring of AP^ 
AP, 
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Proof. (Al) is clear. The involution •* used in (A2) is defined by 



ny) = fiy), yGK, 

for / E AP. (A3) (II) and (A3) (12) follow from the definition of w. (A3) (13) 
follows for example from |16l Theorem 2.6 and Example 2.10], where it is 
shown that w is a topological index on AP, and hence in particular, it is 
locally constant. 

Finally, (A4) follows from [71 Theorem 1, p. 776] which says that / E AP~^ 
satisfies 

inf 1/(^)1 >0 (5.2) 

lm(s)>0 

iff inf \f{y)\ > and w{f) = 0. But 

inf |/(y)| >0 



is equivalent to / being an invertible element of AP by the corona theo- 
rem for AP (see for example Exercise 18, p.24]). Also the equivalence 
of (j5.2p with that of the invertibility of / as an element of AP'^ follows 
from the Arens-Singer corona theorem for AP^ (see for example [3l Theo- 
rems 3.1, 4.3]). Finally, the invertibility of / E i? in i? is equivalent to the 
invertibility of / as an element of AP'^ since is a full subring of AP"^ . □ 

Remark 5.5. Specific examples of such R are AP'^ and APW^. More 
generally, let S C [0, -|-oo) be an additive semigroup (if A,// E S, then 
A + /i E S) and suppose E S. Denote 

APe = {/ E I a(/) C S} 
APW^ = {/ G APW I a{f) C S}. 

Then APj] (respectively APWj]) is a unital Banach subalgebra of AP'^ 
(respectively APW~^). Let Is denote the set of all maps 9 : T, ^ [0,+oo] 
such that 61(0) = and 6'(A + n) = 0{X) + 6'(/i) for all A,;U E S. Examples 
of such maps 9 are the following. If y E [0,+oo), then 9y, defined by 
^yW — A E H, belongs to Is- Another example is 9^, defined as 
follows: 

ifA = 0, 
+00 ifA/0. 

So in this way we can consider [0, -|-oo] as a subset of Is- 

The results [3 Proposition 4.2, Theorem 4.3] say that if Is C [0, +oo], 
and / belongs to Ai-*s (respectively to APWt,), then / belongs to inv APj: 
(respectively to inv APW^;) iff (|5.2p holds. So in this case AP^ and APWy; 
are unital full subalgebras of AP^ . 

In the definition of the i/- metric given in Definition 12.81 corresponding to 
Lemma 15.41 the || • ||oo now means the following: if F E {AP)^^'^, then 

||F||oo = sup|F(2/)|. 



.(A) 
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This follows from (12. 6p , since M is dense in the maximal ideal space DJl (which 
is the Bohr compactification Mb of M) of the Banach algebra S = AP; see 
[m Exercise 18, p. 24]. 



Remark 5.6. It was shown in [2] that AP^ and APW~^ are projective 
free rings. Thus \i R = AP+ or APW"^ , then the set S{R,p,'m) of plants 
possessing a left and a right coprime factorization coincides with the class 
of plants that are stabilizable by [HI Theorem 6.3]. 

5.3. Algebras of Laplace transforms of measures without a singular 
nonatomic part. Let C+ := {s G C | Re(s) > 0} and let denote the 
Banach algebra 

fa€L\0,oo), {fk)k>oe^\ \ 



A+ = \ s{e C+) ^ lis) + fke-''' 



— to < ti,t2,t3, . . . 



K k=0 
equipped with pointwise operations and the norm: 

00 

l|i^ll = ll/aLi + ll(/fc)fc>o||,i, F(s) = /,(s) + ^/fce-^*'= (5GC_ 

fc=0 

Here fa denotes the Laplace transform of fa, given by 

fOO 



poo 

fa{s)= / e"''fa{t)dt, SGCh 

Jo 



Similarly, define the Banach algebra A as follows ([13]): 



00 



fa^L^iR), {fk)kezee\ 

■ ■ ■ ,t-2,t-i <0 = to< ti,t2, ■ ■ ■ 



A=liy{eiR)^ fa{iy) + Y^ fue 

K k=—oo 

equipped with pointwise operations and the norm: 

00 

l|i"ll = ll/alLi + ll(/fc)fcgz||,i, F{iy):=fa{iy)+ ^ 

fc=— 00 

/oo 
e~''^^ fa{t)dt, {y G M) 
-00 

It can be shown that (M) is an ideal of A. 



00 



For F = /, + Ae"* *^ G A, we set Fap(%) = ^ ^ 

fc=— 00 fc=— 00 

11 F = fa + Fap G iiiv A, then it can be shown that FAp{i ) G inv AP as 
follows. First of all, the maximal ideal space of A contains a copy of the 
maximal ideal space of APW in the following manner: if G M[APW), 
then the map ^ : A ^ C defined by = ^{fa + Fap) = ^{FAp{i-)), 

(F G A), belongs to M{A). So if F is invertible in A, in particular for 
every <I> of the type describe above, / ^{F) = Lp{FAp{i-))- Thus by the 
elementary theory of Banach algebras, FAp{i-) is an invertible element of 
AP. 

Moreover, since L^(]R) is an ideal in A, F^pfa is the Fourier transform of 
a function in L^(]R), and so the map y i-> 1 + {FAp{iy))"^fa{iy) = PApiy) 
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has a well-defined winding number w around 0. Define W : inv ^ — )• R x Z 
by 

W{F) = {w{FAp)M^ + FXlfa)), (5.3) 
where F = fa + Fap G inv A, and 

w{Fap) := jim^ ^ arg (F^p(ii?)) - arg {FAp{-iR))^ , 

w(l + F^j,/,):=— (^arg(l + (F^p(fy)) Va(^?/)) j- 



Lemma 5.7. F = fa + Fap £ A is invertible iff for all ?/ G M, F{iy) / 
and inf \FAp{iy)\ > . 



Proof. The 'only if part is clear. We simply show the 'if part below. 
Let F = fa + Fap G ^ be such that F{iy) 7^ for ah ?/ G R and 

inf \FAp{iy)\ > 0. 



Thus FAp{i-) is invertible as an element of ^P. Hence F = FAp{l + faFXp) 
and so it follows that (1 + faF^p)(iy) 7^ for all y G R. But by the corona 
theorem for 

W := LMR) + C 

(see [121 Corollary 1, p. 109]), it follows that 1 + fa^AP invertible as an 
element of W and in particular, also as an element of A. This completes the 
proof. □ 



Lemma 5.8. Let 



R 
S 
G 
i 



a unital full subring of A'^ , 
A, 

R X Z, 
W. 



Then (Al)-(A4) are satisfied. 

Proof. (Al) is clear. The involution •* in (A2) is defined by 



F*{iy) = F{iy), y G R, 

for F & A. (A3) (12) is now easy to see from the definition of W. Also, 
(A3) (II) follows from the definition of W as follows. Let F = fa + Fap and 
G = ga + Gap ■ Then we have 

w{FapGap) = w{Fap) + w{Gap) 
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from the definition of w. Tlius 

W{FG) = W{{fa + FAp){ga + GAp)) 

= WifaQa + TaGAP + QaFAP + FapGap) 

= (w(l + {FApGAPy\Iaga + faGAP + faFAp), w{FapGap)) 

= {^{{l+FX^fa){l + G2l>9a)),w{FAp)+w{GAp)) 
= (w(l + F^^fa) + w(l + G^'pga),w{FAp) + w{Gap)) 

= W{fa + FAP) + W{ga + GAP). 

So (A3) (12) liolds. 

The local constancy of W demanded in (A3) (13) can be seen in the fol- 
lowing manner. We have already noted that w is locally constant on inv AP 
and w is locally constant on inv C(T). Note that ^il + F^pfa) defined above 
is just v{(p) where 

^ 1 -I- e*^ 

^{6) = (1 + Fx'pfa){iy), where iy = j—r^, G (0, 27r). 

Hence (A3) (13) follows. 

Finally we check that (A4) holds. Suppose that F = fa + Fap belonging 
to ^"'"n(inv A), is such that W{F) = 0. Since F is invertible in A, it follows 
that FAp{i-) is invertible as an element of AP. But w{Fap) = 0, and so 
Fap{}-) G AP^ is invertible as an element of AP'^ . But this implies that 
1 + F^p fa belongs to the Banach algebra 

W+ := LHO^)+C. 
Moreover, it is bounded away from on iM since 

1 + ^Ip7a = — , 
^AP 

and F is bounded away from zero on im. Moreover w(l + F^p/a) = 0, 
and so it follows that 1 + F^p fa is invertible as an element of W^, and in 
particular in A'^. Since F = (1 + FXpfa)FAp and we have shown that both 
(1 + F^pfa) as well as Fap are invertible as elements of A'^, it follows that 
F is invertible in A~^. □ 

An example of such a R (besides A~^) is the algebra 

L^'(0^oo)+APWj,{i-) := {fa + FAP : fa G 1^0,+^), FAp{i-) G APWj,}, 



where S is as described in Remark 15.51 

In the definition of the i/- metric given in Definition 12.81 corresponding to 
Lemma 15.81 the || • ||oo now means the following: if F G A^^"^, then 

||F||oo = sup \F{iy)\. 
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This follows from ()2.6p . since M is dense in the maximal ideal space dJl of 
the Banach algebra S = A; see [Ml Theorems 4.20.1 and 4.20.4]. 

Remark 5.9. It was shown in [2] that A'^ is a projective free ring. Thus the 
set §{A'^ ,p,m) of plants possessing a left and a right coprime factorization 
coincides with the class of plants that are stabilizable by jl8| Theorem 6.3]. 

5.4. The polydisk algebra. Let 

B*" := {(zi,... G C" : < 1 for i = 1,... ,n}, 

B" := {(zi,... ,z„) G C" : < 1 for i = 1,... ,n}, 

T" := {(zi,...,z,)EC":|zi| = lfori = l,...,n}. 

The polydisk algebra ^(B") is the set of all functions / : B" — )• C such that 
/ is holomorphic in B" and continuous on B". 

If / G A(B'^), then the function fa defined hy z ^ f{z, . . . , z) : B ^ C 
belongs to the disk algebra j4(D), and in particular also to C(T). The map 

/ ^ Uh^Jd) : ^(O") ^ C7(T") X C(T) 

is a ring homomorphism. This map is also injective, and this is an immediate 
consequence of Cauchy's formula; see [191 P-4-5]. We recall the following 
result; see [El Theorem 4.7.2, p. 87]. 

Proposition 5.10. Suppose that ^ = (V'l, • • • , ^n) is a continuous map 
from B into B", which carries T into T" and the winding number of each 
ipi is positive. Then for every f £ ^(B'"), /(^(B) U T'') = f{W). 

Lemma 5.11. Let 

R = a unital full subring of A(B"'), 
S := C(T") X C(T), 
G := Z, 
L := ^ w(/i)). 

Then (A1)-(A4) are satisfied. 

Proof. (Al) is clear. The involution •* in (A2) is defined as follows: if 

{f,g) G C(T") X C(T), then (f^g)* := {f\g*), where 

/*(zi, . . . ,Zn) = /(zi, . . . ,z„), (zi,...,z„) G T", 
g*{z) = W), ^eT. 

(A3) was proved earlier in Subsection 15.11 Finally, we will show below that 
(A4) holds, following [10]. 

Suppose that / G ^(B") is such that f\r^ G inv C(T"), fa G inv C(T) and 
that w(/rf) = 0. We use Proposition 15.101 with ^{z) := (z, . . . ,z) (z G B). 
Then we know that / will have no zeros in B" if /(^'(B)) does not contain 
0. But since fd G inv C(T) and w(/d) = 0, it follows that fa is invertible as 
an element of ^(B) by the result in Subsection 15.11 But this implies that 
/(^(B)) does not contain 0. 
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Now suppose that / E A(D'') with f\jn G inv C(T"), fd G inv C(T), and 
that it is invertible as an element of A{W"'). But then in particular, fa is an 
invertible element of ^(B), and so again by the result in Subsection 15.11 it 
follows that w(/rf) =0. □ 



Besides y4(B") itself, another example of such an R is i?/f°°(B"), the set of 
all rational functions without poles in B'^. 

In the definition of the z^-metric given in Definition 12.81 corresponding 
to Lemma 15.111 the || • ||oo now means the following: ii F = {G,H) G 
(C(T") X C(T))P^'", then 



||F||oo = max < max \ G{z) \ , max \ H{w) \ >. 

This follows from ()2.6p . since the maximal ideal space 93T of the Banach 
algebra 5 = C(T") x C(T) can be identified with T'' U T. 

Remark 5.12. By [2], it follows that j4(D") is a projective free ring, since 
its maximal ideal space the polydisk B" is contractible. Thus the set 
S(^(B"),p, m) of plants possessing a left and a right coprime factorization 
coincides with the class of plants that are stabilizable by |18l Theorem 6.3]. 

Remark 5.13. Roughly, the index function l: inv 5" — )• G in all the exam- 
ples given above (Sections l5.lH5.4p can be viewed as generalizations of the 
winding number for a continuous nonvanishing function on the unit circle. 
Another important application of such index functions, apart from robust 
control theory as presented here, is to the Fredholm theory of various classes 
of operators (e.g., Toeplitz, Wiener-Hopf, convolution) associated with the 
function. In this context we mention that Murphy [16] has given an ab- 
stract quantized C*-algebra setting which, among other things, unifies the 
connection between analytic index and Predholm index for the C(T)-setting 
of Section 15.11 and the ^P-setting of Section 15. 2[ There has also been a 
substantial amount of other work (see the books [SJ Hj ) where the analytic 
index has been extended to more general classes of functions (e.g. piecewise- 
continuous) in order to develop the Fredholm theory for more general classes 
of Toeplitz operators. On the other hand, the index theory for semi-almost 
periodic symbols (a version of the Callier-Desoer class where fa is only re- 
quired to be continuous on the extended imaginary line and where f — fa 
is required only to be AP rather than APW) follows a different more com- 
plicated path rather than making use of the index function W as in ()5.3p . 
Similarly, the Fredholm theory for Toeplitz operators on the quarter plane 
(associated with continuous functions on the bitorus T^) (see [6l Chapter 8]) 
makes use of the Z^-valued index associated with the winding number of a 
function / on taken with respect to each variable separately, rather than 
with the index l as in Lemma l5.11[ 
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6. Further directions 

It was shown in |25j that when R comprised rational functions without poles 
in the closed unit disk, then the bound established in Theorem 14.61 is the 
best possible one in the following sense: 

(P'): C satisfying ^Pa,c > stabilizes P only if du{P, Pq) < m. 

Since this property of already holds in the rational case, we expect the 
same to hold also in the specific examples considered in the previous section. 
We leave the question of investigation of whether the property (P') always 
holds in our abstract setup for future work. 
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